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Abstract We searched for X-ray bursts in XMM-Newton archival data of X-ray sources in M 31 globular clusters (GCs) and 
GC candidates. We detected two bursts simultaneously in EPIC pn and MOS detectors and some more candidates in EPIC pn. 
The energy distribution of the burst photons and the intrinsic luminosity during the peak of the bursts indicate that at least the 
strongest burst was a type I radius expansion burst. The bursts identify the sources as neutron star low mass X-ray binaries in 
M 3 1 . The type I X-ray bursts in M 3 1 are the first detected outside the Milky Way and show that with the help of XMM-Newton 
X-ray bursts can be used to classify neutron star low mass X-ray binaries in Local Group galaxies. 
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1. Introduction 

Type I X-ray bursts (hereafter X-ray bursts) were detected in 
low mass X-ray binaries (LMXBs) in the Milky Way and are 
explained as thermonuclear runaways in freshly accreted mate- 
rial on the surface of a neutron star They may re-occur on a time 
scale of hours to days, when enough new material is accreted, 
compressed and heated, to densities and temperatu res adequate 
for another thermonuclear ignition ( for reviews see lLewin et alJ 
ll993VStrohmaver & Bildsten'2003V X-ray bursts show a large 
variety in profiles and peak flux. The persistent luminosities of 
the sources when bursting, vary from ~0.2 to well bellow 0.01 
of the Eddington luminosity. Burst rise times last from less than 
a second to ~10 s, and decay times are in the range of ~10 s 
to minutes. In general, burst profiles depend strongly on pho- 
ton energy, with decays which are much shorter at high photon 
energies than at low energies. This energy dependence of the 
burst profile corresponds to a softening of the burst spectrum 
during decay, which is the result of cooling of the neutron star 
photosphere. The time dependent energy spectra can be mod- 
eled by black body radiation of changing temperature. 

During some bursts the energy release is high enough that 
the luminosity at the surface of the neutron star reaches the 
Eddington limit leading to an expansion of the neutron star at- 
mosphere (radius expansion X-ray bursts). These bursts show 
an energy dependent double peak structure. During expansion 
and subsequent contraction the luminosity remains almost con- 
stant near the Eddington limit. However, the black body tem- 
perature may quickly drop from more than 3 keV to well below 
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0.5 keV when the radius of the photosphere expands before 
it slowly rises again during contraction to 2-3 keV and then 
cools thereafter with continuously decreasing luminosity. For 
a detailed description of the temperature and flux development 
of several radius expansion bursts from 4U/MXB 1820-30 see 
Haberl et al. ( 1987). The radius expansion bursts of about two 
thirds of the sources reach a critical bolometric luminosity of 
~3.8xlO^^ erg s"', corresponding to the Eddington limit for 
hydrogen-poor matter and may be used as standard candles 
(see e.g. lKuulkers et al. 2003). The expansion phase may last 
for up to 30 s. To date thirteen luminous globular cluster (GC) 
X-ray sources are known in the Milky Way, twelve of which 
show X-ray bursts (and sometimes even radius expansion X-ray 
bursts) and therefore are believed to be LMXBs containing a 
neutron star (see e.g. Verbunt & Lewin 2004). The distances of 
the sources vary between about 4 and 12 kpc. Up to now, no 
X-ray bursts have been detected outside the Milky Way. 

With the high collecting power and good position reso- 
lution of the XMM-Newton EPIC cameras we now have the 
possibility to detect X-ray burst sources in all galaxies of the 
Local Group. The X-ray burst signature then will identify these 
sources as neutron star LMXBs. Using the maximum burst flux 
as standard candle would allow to determine the distance of 
the burst source. However, due to the limiting statistics, dis- 
tance determinations by this method can not compete with 
the accuracy of other methods. The persistent luminosity of 
many bursting LMXBs may be below the limiting point source 
sensitivity of current observations. Therefore, to find all X- 
ray burst sources in neighboring galaxies, one would have to 
search the entire extent of the galaxy with the resolution of 
the instrument point spread function (PSF). A less complete 
method would be a search in LMXB candidates selected by 
other means (e.g. as GC X-ray sources). Radius expansion 
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Table 1. XMM-Newton EPIC count rates for the pn and one 
MOS instrument (thin filter full frame mode, 15 arcsec on-axis 
extraction radius corresponding to approximately the half en- 
ergy width of the PSF) in different energy bands (in keV) ex- 
pected from a source emitting a black body spectrum with tem- 
perature Tbb and intrinsic bolometric luminosity of lO"'^ erg s"' 
at the M 3 1 distance assuming Galactic foreground absorption 
of 7x1020 cm-2. 



EPIC pn EPIC MOS 

0.2-0.5 0.5-4.5 4.5-7.0 0.2-0.5 0.5^.5 4.5-7.0 

(keV) (IQ-^ counts s"') 
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124 


125 
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30 





0.2 


70 
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333 
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bursts would be most easily picked up as they are the most 
luminous "standard" type I X-ray bursts. Besides these bursts 
one may detect type I "superbursts" lasting for several hours 
with peak luminosities above lO-'^ erg s"' and maximum black 
body temperatures of 2-3 keV similar to the ones detected with 
Be ppoSAX and RXTE from LM XBs in the Milky Way (see 
e.g.'Stro hmaver & Bildstenl2003l) . 

An ideal target for a search for X-ray burst sourc es is the 
brigh t Local Group spiral M 31 (distance 780 kpc, iHollanj 
1 1998 ': 'Stane k & Garnavichl[l998! ) with its moderate Galactic 
foreground absorption (A^h= 7x10-" cm"^, .Stark et al..,1992 ). 
Many GC a r id GC candidates are known in M 31 (see e.g. 
iGalleti et al." '2004') which may host LMXBs. IPietsch et alJ 
t2005., hereafter PFH2005) prepared a catalogue of M 3 1 point- 
like X-ray sources analyzing all observations in the XMM- 
Newton archive which overlap at least in part with the optical 
D25 extent of the galaxy. In total, they detected 856 sources and 
identified and classified them using their X-ray properties and 
correlations with sources in other wavelength regimes. Twenty- 
seven of the X-ray sources are identified with GCs and ten 
with GC candidates. They are - extrapolating from bright GC 
sources in the Milky Way - most likely neutron star LMXBs. 
In this paper we report on a search for X-ray bursts in the M 3 1 
XMM-Newton data of these LMXB candidates. 

2. XMM-Newton observations and analysis 

We searched f or M 31 bu rst sources in the sam e archival 
XMM-Newton djansen et al.. 2001 ) EPIC dStriider et al. 2001: 
iTurner et alJ200l|) observations that were used by PFH2005 for 
the creation of the M 31 source catalogue, i.e. pointings cl to 
c4 to the galaxy centre, nl to n3 to the northern disk, si and s2 
to the southern disk and h4 to the northwest halo (see Table 1 of 
PFH2005 for details). The observations were performed in the 
fuU frame mode (time resolution for pn and MOS 73.4 ms and 



2.6 s, respectively) using medium or thin filter with low back- 
ground exposure times of about 10 to 50 ks. For the light curve 
analysis we used event pattern up to 12 (singles to quadruples). 
We restricted our search to the stringent low background times 
used by PFH2005 to avoid spurious burst candidates due to 
background flares which would show up in all EPIC instru- 
ments at the same time and might be misidentified as bursts. 

To get an estimate for the counts that can be expected in 
the EPIC detectors from an X-ray burst in M 31, we calcu- 
lated count rates for typical burst spectra (black body spectra 
with temperatures of 0. 1 to 3 keV) for a source with an intrin- 
sic luminosity of lO''^ erg s"' (Table 0. It is clear that only 
bursts with radius expansion, that radiate at the Eddington lu- 
minosity for at least several seconds will give enough photons 
(more than about 10) to be significantly detected in pn detector 
and possibly also several photons for burst confirmation in the 
MOS detectors. The expected number of photons in the indi- 
vidual MOS instruments is lower than for the pn by a factor of 
two or more depending on the assumed spectrum (see TableQ. 
The drop of the black body temperature during the expansion 
to values below 1 keV additionally favors their detection as it 
shifts the main energy output into the band where EPIC is most 
sensitive. Taking all the points mentioned above into account, a 
search strategy for bursts should be based on pn light curves of 
candidate sources with time resolutions of about 10 s or longer 
and use MOS count rates for confirmation. 

Therefore, we created light curves for EPIC pn in the 0.5- 
4.5 keV band with 500 s, 100 s and 20 s resolution for the 37 
sources associated with a GC or possibly associated with a GC 
in the catalogue of PFH2005 selecting all events within a radius 
around the source corresponding to the half energy width of 
the PSF at the off'-axis angle. We rejected light curves when the 
source was outside the EPIC pn field of view (FOV) or partly 
covered by EPIC pn CCD gaps. The effective integrated expo- 
sure time for all sources adds up to 2.88x10^ s. We inspected 
the light curves by eye for burst candidates. Two bursts are con- 
firmed in MOS light curves and are therefore regarded as se- 
cure detections. Background flares are ruled out as they would 
appear simultaneously in other source light curves. The bursts 
most likely showed radius expansion and will be discussed in 
greater detail in the following sections. Several other pn burst 
candidates were not in the FOV of the MOS detectors or not 
detected by them. No burst with characteristics similar to the 
superbursts in the Galaxy could be detected. However, as also 
in the Galaxy superbursts are very rare events, it is not surpris- 
ing that none appears in the available M 31 observations. 

The data analysis was performed using tools in the 
XMM-Newton Science Analysis System (SAS) v6.0.0, 
EXSAS/MIDAS 03OCTXXP and FTOOLS v5.2 software 
packages, the imaging application DS9 v3.0b6 together 
with the funtools package, the mission count rate simulator 
WebPIMMS v3.6a and the spectral analysis software XSPEC 
vll.3.1. 



M31 burst source [PFH2005] 253 
XMM-Newton, EPIC 
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Figure 1. XMM-Newton EPIC light curves of source 
[PFH2005] 253 during the M 31 centre observation c4 (OBSID 
01 12570101) on January 6/7, 2002 integrated over 120 s. Time 
zero corresponds to the start of the observation. 



M31 burst source [PFH2005] 253 



XMM-Newton. EPIC 
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Figure!. Combined XMM-Newton EPIC Ught curve (pn, 
MOSl and MOS2 added) of the burst in Fig [T] of source 
[PFH2005] 253 integrated over 10 s. The dashed histogram 
shows the light curve of the 0.2-4.5 keV events. 



3. The X-ray burst source in the globular cluster 
[WSB85] S5 15 

The brightest M 31 burst in the EPIC pn and both MOS in- 
struments was detected from source [PFH2005] 253 in GC 
[WSB85] S5 15 (.Wirthet al. 19 85) during observation c4 start- 
ing at UT 6: 14:50 on January 7, 2002 (Fig.[T]and|3. It reached 
its maximum luminosity within 10 s, stayed at the maximum 
(count rate in 0.2-7 keV band of (1.4+0.4) ct s"') for less than 
20 s and decayed to background level within about 150 s. Due 
to limiting statistics the structure of the bursts can not be fur- 
ther resolved. In total, more than 50 burst counts are detected. 
As can be seen from Fig.|2lphotons in the 4.5-7 keV band are 
only detected after the maximum of the burst as expected for 
the temperature development in a radius expansion burst. Also 
the maximum count rate corresponds to that of a M 3 1 source 
with a 1 keV black body spectrum radiating at 3.8x10^^ erg s"', 
i.e. at the Eddington limit for hydrogen-poor matter. 




1 2 
Channel Energy (keV) 

Figures. XMM-Newton EPIC spectrum of source [PFH2005] 
253 for observation c4. The best fit to the three EPIC spectra 
- absorbed bremsstrahlung plus black body model - is shown 
in the upper panel (pn is the upper spectrum and the two MOS 
spectra overlap). 



Source [PFH2005] 253 was in the EPIC FOV during the 
four XMM-Newton M 31 centre pointings cl to c4 (average 
absorbed luminosity of 3.4x10-^^ erg s"' in 0.2-4.5 keV band). 
Between the pointings which are equally spread over two years, 
its average luminosity varied by a factor of ~6. During pointing 
cl [PFH2005] 253 was partly hidden in pn by a CCD gap. We 
therefore only used pointings c2, c3 and c4 for the burst search 
(effective low background time of 75.7 ks). 

During observation c4, more than 2 500 photons were 
detected in the EPIC detectors from [PFH2005] 253. An 
absorbed bremsstrahlung plus black body model describes 
the persistent energy spectrum (see Fig. |3j best fitting 
model pai-ameters: A^h= (T_l^)xlO^^ cm^^, Tbrems = 5+?° keV, 



<rems ~ _2 

Tbb = 0.18!^;^^ keV, reduced ;tL = 1-66 for 45 degrees of free- 
dom) and corresponds to an absorbed luminosity and intrin- 
sic luminosity of the source in the 0.2^.5 keV band of 4.1 
and 5.6x10^^ ergs respectively. The source was also de- 
tected by Chandra ACIS-I and HRC as a variable source (r3- 



41. Kong et al.. (.2002.) : CXOM31 J004221.5H-41 1419, 
ll2002h l. 



rce (r3- 
.Kaaretl 



4. The X-ray burst source in the globular cluster 
B150 

The EPIC pn and the MOS2 instruments detected a burst from 
M 31 source [PFH2005] 396 in GC B150 (G alleti et al. 2004) 
during observation c3 starting at UT 16:58:00 on June 29, 2001 
(Fig. 13. In MOSl the source was heavily obscured by a CCD 
gap and therefore MOS 1 data could not be used for the analy- 
sis. The burst was significantly fainter (about a factor of three) 
than the burst of [PFH2005] 253 and only about 15 counts are 
detected, four or them in M0S2. The burst lasted for about 40 s 
with a rise to maximum within 10 s. Due to limiting statistics 
the structure of the bursts can not be further resolved. The count 
rate at maximum provides us with a lower limit to the peak lu- 



4 



W. Pietsch and F. Haberl: Type I X-ray bursts in M 3 1 



M31 burst source [PFHZOOo] 396 
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Figure 4. Source [PFH2005] 396: XMM-Newton EPIC pn and 
MOS2 light curves during the M 31 centre observation c3 
(OBSID 0112570701) on June 29, 2001 integrated over 120 s 
(left). Combined XMM-Newton EPIC light curve (pn and 
MOS2 added) of the burst integrated over 10 s (right). Time 
zero corresponds to the start of the observation. 



for some time during the contraction phase of the radius ex- 
pansion when the burst is radiating at the Eddington luminos- 
ity. Of specific importance is the possibility to simultaneously 
detect a burst in more than one EPIC instrument, which gives 
high confidence in the credibility of a burst even if only few 
photons are involved. In this work only known LMXB candi- 
dates in M 3 1 from PFH2005 were investigated. The success 
encourages a more systematic survey covering the total area of 
M 3 1 and also other Local Group galaxies to identify neutron 
star LMXBs just from their bursting properties. However, to 
be able to use the bu rsts as standard candles as proposed by 
iKuulkers et all ( l2003t) even better statistics at higher time reso- 
lution are needed to analyze the burst signature in more detail. 
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minosity in M 3 1 (assuming a 1 keV black body spectrum) of 
about 10-"* erg s"'. 

Source [PFH2005] 396 was in the EPIC FOV during the 
four XMM-Newton M 31 centre pointings cl to c4 (average 
absorbed luminosity of 1.7x10-^^ erg s"' in 0.2^.5 keV band). 
Between the pointings, its absorbed luminosity varied by a fac- 
tor of ~1.7. We used pointings cl, c2, c3 and c4 for the burst 
search (effective low background time of 99.3 ks). 

During observation c3, the source was at a luminosity of 
(2.0+0.2)xl0^^ ergs"^ The faintness of the source and the 
low background exposure time for observation c3 that was a 
factor of two shorter than for observation c4, prevented a sim- 
ilar spectral modeling as for [PFH2005] 253. The source was 
also detected by Chandra ACIS-I and HRC as a variable source 
(r3-18 Ikongetal. L2002) : CXOM31 J004307.4-H4 12021, 
lKaareill2002h l. 

5. Discussion 

We have detected the first type I X-ray bursts outside the 
Milky Way from two M 31 GC source candidates of PFH2005 
which showed persistent luminosities while bursting of below 
~5xlO^^ erg s"'. In search for bursts we checked low back- 
ground EPIC pn light curves (covering in total 2.88x10^ s) of 
sources associated with GCs or possibly associated with GCs. 
The persistent luminosities of part of the sources was above 0.2 
times the Eddington luminosity of a neutron star and therefore, 
extrapolating from the experience of Galactic LMXBs, were 
not expected to show bursts during these bright states. No burst 
with characteristics similar to the superbursts in the Galaxy was 
detected. 

The peak luminosities of the two bursts in [WSB85] S5 15 
and B150 are comparable to those of bright radius expansion 
bursts in the Milky Way. The detections demonstrate the pos- 
sibility to classify neutron star LMXBs in the Local Group 
by their bursting signature using just X-ray information from 
XMM-Newton EPIC. The energy response of the EPIC instru- 
ments is most sensitive to detect radiation with black body tem- 
peratures in the range 0.2-1 keV which are normally covered 
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